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Abstract 

We study the influence of short range antiferromagnetic correlations between local 
/-electrons on the transport and thermodynamic properties of Kondo insulators. We 
use a Kondo lattice model with an additional Heisenberg interaction between nearest 
neighboring f-moments, as first proposed by Coqblin et al.[l] for metallic heavy 
fermions. The inter-site magnetic correlations produce an effective bandwidth for 
the /-electrons. They are treated on the same footing as the local Kondo correlations 
such that two energy scales appear in our approach. We discuss the competition 
between these two scales on the physical properties. 
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1 Introduction 



Kondo insulators form a group of compounds that at high temperatures behave 
as dirty metals but at very low temperatures have their thermodynamic and 
transport properties determined by the existence of a small gap (10-lOOK) that 
arises from the hybridization of local electrons and the conduction band [2,3]. 
This family of compounds can be characterized as strongly correlated semi- 
conductors due to the f or d character of the relevant electrons and includes 
FeSi [4], Ce^BuPh [5-7], SmB^ [8], YhBi2 [9] and CeFe^Pu [10]. 

Many theoretical models have been used to describe Kondo insulators. Some 
of them considered a two-band system, a large uncorrelated band of s electrons 
and another narrow correlated band which describes the f or d electrons [11- 
13]. Others considered an Anderson lattice, a conduction band and localized / 
states, with the correlation in these states generally treated within the slave- 
boson method [14,15]. 
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Although Kondo insulators do not present long range magnetic order, the oc- 
currence of short range antifcrromagnetic correlations in these materials has 
been observed by inelastic neutron scattering [16], suggesting that a good 
model to describe Kondo insulators must include antiferromagnetic correla- 
tions between neighbors /-moments. We also point out that the magnetic 
susceptibility of Kondo insulators in the Curie- Weiss regime has a negative 
Curie temperature, as shown in Table 1 for some compounds. This is a clear 
indication of the presence of antiferromagnetic correlations in these materials. 

In this letter we study some physical properties of the Kondo lattice taking 
into account short-range antiferromagnetic correlations between the localized 
/-electrons (only first neighbors) as first proposed by Coqblin et aL[l]. These 
authors have used this approach to describe metallic heavy-fermions while 
here we shall use it to investigate Kondo insulators. We calculate the density 
of states and study its behavior with increasing temperature and increasing 
magnetic correlations. We evaluate some physical properties such as optical 
conductivity, magnetic susceptibility and electrical resistivity. The format of 
the present letter is as follows. In section 2 we present the Hamiltonian and 
describe the approach used. In section 3 we calculate the density of states 
and some physical quantities in order to analyze the influence of magnetic 
correlations in the system. 



2 Model Hamiltonian 

We consider the following Hamiltonian to describe the system [1] : 

H ^knl^ + ^0 X! ^i<y - JK^Si - Si - Jh ■ 6^+5, (1) 

k,a i,(T i i,5 

where n^^ = c\^Cka, n{a = fLficr- ^^e spin operators associated to 

localized /-moments and si to conduction electrons at site i. The flrst term 
represents the conduction band. The operators c|.^ and Cka, respectively, create 
and annihilate electrons in the conduction band state labeled by the wave vec- 
tor k and spin a. The second term, proportional to Eq, represents the binding 



Compound 


e 


Ce3Bi4Pt3 

YbBu 

FeSi 


- 125 K [7] 
-79 K [17] 
-1030 K [4] 



Table 1 

Table of Curie Temperature for some compounds 
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energy of the /-electrons. The third term is the s — / exchange interaction Jk 
which gives rise to the Kondo couphng. In this model, following Coqblin et 
al. [1], we add a Heisenberg-like interaction Jh between nearest neighboring 
/-moments. 

We examine the model above, specifically in the situation where the ground 
state of the system is an insulating state. First we construct an effective Hamil- 
tonian, rewriting the original one, Eq.(l), in terms of new operators, as per- 
formed by Ruppenthal et. al [18]. Since we intend to describe the Kondo 
effect, we choose an operator Aj^o- that couples c and / electrons. In order to 
describe short-range magnetic correlations we introduce another operator Fj 
that couples / electrons in neighboring sites. These Hermitian operators are 
given by: 

= ~7^('-'icr/io- + flakier) (2) 



and 



I 

71 



{fitjfi+SjO- + fi+S,afi' 



(3) 



The spin components can be written in terms of the new operators. Using an 
additional constraint, that excludes the double /-occupancy, we get: 



(4) 



- Vf2 



where nf — nf^ + nf^. Since we are describing an insulating state at low tem- 
peratures, we use the constraints < n{ >—< >= 1. Furthermore applying 
the decoupling, 

< rCn% >^< < > <+ <n%>nl, (5) 



we find that the terms containing the number operators cancel out. Note that 
we are interested in non-magnetic solutions such that < nf^ >=< nf_^ >. 
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Finally we can write the Kondo and Heisenberg parts of the Hamiltonian as: 

Hk — -^Jr X] O^iT + K,-a)K(n (6) 
ia 

Hh = -Jh X! (^1+5,(7 + ^i+5-a)^i+5,a- (7) 
i+6(T 



Consistently with the decoupling used above, Eq.(5), we deal with the product 

of operators in equations (6) and (7) introducing the mean fields A, =< Aj^o- > 
and Vi =< ti^s,a >■ Furthermore due to translational invariance, Aj = A and 
Fi = r. Performing a Fourier transform, using £k = (kR) where z 

is the number of neighbors and 2W is the band width, we obtain the effective 
Hamiltonian: 



H' = Y: ekanl, + JkX E i4afka + fUka) + E (^0 - te/c)nL (8) 

k,a kcr k,cr 

— QiJj^y? — 2 JffF^, 

where 

6 = ^r. (9) 

We have then arrived at a new picture where the magnetic correlations Jh give 
rise to a narrow band of /-electrons of effective width 261^ (b is temperature 
dependent). The interaction Jk (using Jr- < 0, since we are describing the 
Kondo effect) introduces an effective hybridization between the conduction 
and the renormalized / bands. 

The new quasi-particles associated with the Hamiltonian (9) can now be ob- 
tained. This is done by calculating the one-electron Green's functions which 
are found from their equations of motion. These Green's functions are given 

by, 

g!M - - (10) 
G^M = ^^^^^^^f Vt). (11) 
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where =<< al^\Pk'a » ^or a, P — c, , f and 

oj - {bek + Eq) 



UJi{ek) -CU2(£fc)' 



(13) 
(14) 



We have expressed the Green's functions in terms of simple poles correspond- 
ing to two hybridized quasi-particle bands with dispersion relations given by 

1 1 

c^i;2(£fe) = ^[{l + b)ek + £^0] ± -^l[{l-h)sk-Eof + AJl\^. (15) 





J„=0 




J„<0 













Fig. 1. The energy dispersion relations 0Ji-^2{k) for the hybridized bands found in 
the mean-field approximation using = —Wcosika/Ti), where G/2 = t: /2a is half 
a reciprocal lattice vector. Solid line is used for Jh = and dashed line for Jh < 0. 

As shown in figure 1 the two bands are separated by an indirect gap between 
the lower band at A; = G/2 and the upper band at k = 0. For increasing values 
of the magnetic correlation \Jh\, the indirect gap gets larger while the direct 
gap remains constant. The direct gap A^ir is governed only by Jk- It is easy 
to see that A^ij. — 2JkX and it occurs at Eq — (1 — b)ek- 

An insulating ground state, as observed in Ce3Bi4Pts[5] or YbBi2[Q]: is ob- 
tained when the lower band is completely filled. This corresponds to the half- 
filled band case where the occupation numbers assume the values < n{ >=< 
n1 >= 1. Furthermore, we take Eq = /j, = and symmetric bands with 
respect to the chemical potential jj,. Since the /-moments are correlated anti- 
ferromagnetically, which corresponds to < 0, for Jk ^ we always have 
an insulating ground state with the chemical potential fixed in the middle of 
the gap of the density of states. 

The parameters A and F are obtained from the two coupled equations that 
have to be solved self-consistently for each temperature: 
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^ - UJ2) \uj -cui ijj -IJJ2) 

^ - a;2) V a; - a;i a; - a;2 / 

where /(cu) = exp/3(u!-/^)+i Fermi distribution and ci;i;2 = '^\-;2[Bh) is given 

by Eq.(15) 

The present method has many similarities with the slave-boson mean field 
approach and like it, presents a critical temperature where the conduction 
and the / electrons become decoupled [14,15]. In our model, this arises for 
A(T'c) = r(Tc) = which implies the vanishing of the effective hybridization 
term A. So, these mean-field methods are useful to study the properties of 
Kondo insulators only in the low temperature regime, below the critical tem- 
perature Tc of the spurious phase transition. We shall present results here for 
T <C Tc such that the present approach is vahd. 



3 Theoretical Results and Compcirison with Experiments 



As discussed previously, the magnetic correlations modify the band structure 
of the system, increasing the indirect gap. In this section we analyze the in- 
fiuence of these correlations in some physical properties, calculating A and F 
self-consistently for each value of temperature. For simplicity we consider a 
uniform density of states of width 2W for the conduction electrons. 



3.1 Density of states 



Using that < n'^'^ >= -Jm[G^'^(a;)] we obtain the c and / contributions to 
the density of states. In figures 2 and 3, we show the density of states which 
consists of two hybridized bands separated by a gap. The density of states is 
very sharp near the band edges independently of the form of the unperturbed 
bands. 

Temperature renormalizes the density of states decreasing the gap (figure 2). 
With increasing temperatures the density of states become more peaked near 
the band edges. The density of states also depends on the magnetic correlations 
represented here by the parameter b. The last term in Eq.(9) produces an 
effective /-band and its weight renormalizes with b. In a system where the 
Fermi level is inside the gap, the main effects of short range antiferromagnetic 
correlations are to increase the gap and enhance the density of states near the 
band edges, as we can see in figure 3. 
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Fig. 2. / contribution to the density of states for different temperatures. 
Jj^/W = -0.4, Jh/W = -0.1. 
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Fig. 3. / contribution to the density of states for two different values of Jh and 
ksT/W = 0. 

3.2 Optical conductivity 

The optical conductivity can be written as [19]: 

= ^T.I<l-l rfa,pJ(.,a,)p=(.,. + .)liM:i^fci±i^ (18) 

where p^{s, cu) is the one particle spectral density of the conduction electrons. 
In the optical conductivity, shown in figure 4, the behavior of the gap in- 
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Fig. 4. Optical conductivity for Jk/W = —1.0, Jh/W = —0.1 varying the temper- 
ature. The frequency hu is in units of the bandwidth W. 

volves two characteristic temperature scales as observed experimentally for 
Ce^Bi^^Pi-^ [20] and FeSi [4]. Firstly the high temperature regime, compared 
to the width of the conductivity gap Ac, where the gap Ac itself is strongly 
renormalized with T, as we can see in figure 4. The optical conductivity also 
presents a Drude-like feature at low frequencies, as observed recently by Gor- 
shunov et al [21] in dynamical conductivity measurements in StuBq. This peak 
indicates that the low-energy transport is determined by free charge carriers. 
Then the second temperature scale is set by the width of this Drude peak. Only 
at small temperatures the gap is completely opened. At these low tempera- 
tures the variation in the intensity of this peak, whose weight is transferred to 
the high frequency region of the spectrum, is the main effect of temperature. 
The gap itself remaining unchanged as the low-energy behavior changes from 
being dominated by free carriers to more localized carriers. We point out that 
the interaction Jh has no effect in the gap size of the optical conductivity. 



3.3 Magnetic susceptibility and electrical resistivity 

In this subsection we calculate the electrical resistivity and the magnetic sus- 
ceptibility for different values of Jh to analyze the influence of antiferromag- 
netic correlations in these properties. The electric conductivity is obtained 
from the limit u ^ of the optical conductivity. 
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where /3 = 1/T. The resistivity is given by p = 1/(t(0). Notice that this ex- 
pression must be used with care. In the problem considered here, there is 
translational invariance and consequently k is a good quantum number. How- 
ever in real systems impurity scattering is always present and this limits the 
electron mean free path. In the calculations presented below, this is taken 
into account by including a finite lifetime for the conduction electrons. For- 
mally this is done by replacing, uo uo + iF, in the Green 's functions, where 
F is temperature independent. We point out that here we are interested in 
the temperature dependence of the conductivity which is not affected by the 
magnitude of F. 

If an external magnetic field h with its direction along the z-axis is applied, 
the quasi-particle energies become 

i^iA^k) = ^i(£fe) - f^h. (20) 

The free energy in the presence of the field, using the new quasi-particles 
energies (20), is written as 




and the magnetic susceptibility is given by 




Fig. 5. Temperature behavior of the calculated susceptibility (a) and resistivity (b) 
showing activation behavior. The fitting using expressions given in the text are 
shown as full lines for different values of Jh and Jk /W = —0.4. 

The magnetic susceptibility and electrical resistivity can be obtained from the 
expressions above for different values of Jh and Jk fixed. The calculated low 
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temperature resistivity and susceptibility curves were fitted using the activated 
forms x(T) = (C/T) cxp (— A^/ZcT) and p(T) = p^exp {—Ap/kT), wliich de- 
scribe tlie experimental data in Kondo insulators. As we can see in figure 
5, these analytical expressions give also a good description of the theoretical 
results, indicating that our model describes very well the low temperature 
properties of these systems. The fittings of the theoretical results yield val- 
ues for the transport (Ap) and magnetic (A^) gaps [4]. In figure 6 we show 
these gaps for different values of Jh- A^ and Ap increase for increasing values 
of I Jij/Vr| but their ratio remains the same (A^/Ap pa 2.2) suggesting that 
magnetic correlations do not influence the relation between them. 

0.06 



0.04 

< 

0.02 



0.00 

0.00 0.05 0.10 0.15 

Uh/wi 

Fig. 6. Magnetic (full squares) and transport (full circles) gaps versus Jh for 
Jj^/W = -0.4. 

3.4 Analysis 

Experimental results yield different values for the gaps of Kondo insulators. 
The gap measured in optical conductivity Ac has larger values than the gaps 
measured in neutron scattering and susceptibility measurements, A^ or in 
transport measurements Ap[20]. On the other hand, comparing the spin and 
transport gaps, experimental results show that A^ > Ap [4,7]. The present 
theory accounts for these observations and introduces a new discussion about 
the roll of short-range magnetic correlations in the properties of these gaps. 

We have discussed in section 2 the influence of antiferromagnetic correlations 
in the indirect and direct gaps in the density of states. Although in the cal- 
culations of section 3 we used a square band, without a defined dispersion 
relation, we point out that, since our results show that the frequency gap Ac 
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in the optical conductivity does not renormalize as Jh varies, wc may conclude 
that this quantity is determined by the direct gap. This is expected, since due 
to the negligeable momentum of the photon, optical excitations involve essen- 
tially an energy transfer. On the other hand the results in figure (6) show that 
the transport and susceptibility gaps are strongly renormalized by magnetic 
correlations indicating the indirect nature of the gap that determines these 
physical properties. 



4 Conclusions 



In this work we have discussed a model for Kondo insulators which takes 
into account the influence of short-range magnetic correlations in these sys- 
tems. We used a Kondo lattice with an additional Heisenberg term and we 
treated the problem within a mean-field approach in order to obtain an ef- 
fective Hamiltonian where the additional term gives rise to an effective band 
width for the / electrons. Although the Kondo interaction leads to an indirect 
coupling between the local electrons, we have added explicitly to the Hamilto- 
nian an interaction between them. In this way we can make a straightforward 
analysis of the infiuence of these magnetic correlations on the properties of 
the system. 

We have obtained the density of states of the new quasi-particles and ob- 
served that their energy bands are renormalized by temperature and antifer- 
romagnetic correlations. We have calculated the magnetic susceptibility and 
electric resistivity to get the magnetic and transport gaps which turn out 
to be strongly renormalized by magnetic correlations. Calculating the optical 
conductivity we find a Drude-like peak at very low frequencies, as observed 
recently in dynamical conductivity measurements of SmB^^ll]. However, the 
gap in the optical conductivity remains the same for different values of the 
strength of magnetic correlations. This gap is determined by the direct gap in 
the band structure, while the magnetic and transport gaps arise mostly from 
the indirect gap. The antiferromagnetic correlations change substantially the 
physical properties related to the indirect gap but do not change the quanti- 
ties related to the direct one. These observations can produce an insight for 
future experimental analysis. 

This letter discusses explicitly the short-range magnetic correlations which 
have not been explored previously in the theoretical study of Kondo insula- 
tors. The model is able to reproduce their physical properties, including their 
temperature dependence, and is in good qualitative agreement with experi- 
mental results. 
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